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We propose new chaos indicators for systems with extremely small positive Lyapunov exponents. 

These chaos indicators can firstly detect a sharp transition between the Arnold diffusion regime 
and the Chirikov diffusion regime of the Froeschle map and secondly detect chaoticity in systems 
with zero Lyapunov exponent such as the Boole transformation and the S'-unimodal function to 
characterize sub-exponential diffusions. 
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Introduction In weakly chaotic systems with ex¬ 
tremely small Lyapunov exponents, it is well-known that 
it takes a very long time to estimate the largest Lyapunov 
exponent in the order of which is inversely propositional 
of the largest Lyapunov exponent. Thus, there is a prac¬ 
tice that it can take more than ten times longer than the 
Lyapunov time [1], For example, G. J. Sussman and J. 
Wisdom [2] numerically showed the Lyapunov time of the 
solar system is approximately 5 x 10 6 years by comput¬ 
ing over 1 x 10 s years. For investigating nearly integrable 
systems with such weak chaotic property, Froeschle et al. 
proposed a chaos indicator called Fast Lyapunov Indica¬ 
tor (FLI) [ , 3]. If an initial point belongs to a chaotic 
domain, the time evolution of FLI grows linearly. On the 
contrary, if an initial point belongs to a torus domain, the 
time evolution of FLI grows logarithmically [4]. Besides 
the fact that the original key concept of FLI has not been 
changed, OFLI [5] and OFLI 2 [6] are proposed as the 
improvements of FLI, which can reduce the dependency 
of direction of initial variational vectors. In addition to 
nearly integrable systems, in infinite ergodic systems with 
zero Lyapunov exponents, the sub-exponential behavior 
attracts lots of interests [7, 8]. Akimoto et al. [8] pro¬ 
posed generalized Lyapunov exponent which character¬ 
izes super-exponential chaos and sub-exponential chaos. 
That is because the finite time Lyapunov exponent con¬ 
verges in distribution [7, 9]. If an orbital expansion rate 
A grows in order of A ~ exp (n a ), 0 < a < 1, we need to 
get index a to applying the general Lyapunov exponent. 


In this Communication, we propose a new chaos in¬ 
dicator that can detect chaoticity of weak chaotic sys¬ 
tems with extremely small positive Lyapunov exponent 
more rapidly than these existing methods FLI, OFLI 
and OFLI 2 . In addition, this new chaos indicator can 
firstly detect a sharp transition between Arnold diffu¬ 
sion and Chirikov diffusion. Then, we propose another 
new indicator which can find the index a about sub¬ 
exponential chaos whose orbital expansion rate A grows 
A ~ exp(n“), 0 < a < 1 occuring in the Boole transfor¬ 
mation [7, 10] and the S'-unimodal function [11]. 


Ultra Fast Lyapunov Indicator We assume such a dy¬ 
namical system as 
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where H is a Hessian matrix, fi{x n ) is an *th component 
of f(x n ) and e, is a unit vector about ith component. 
Usually, we ignore the terms whose order is greater than 
one and use such a variational equation as 


<Wi = Df(x n )8 n . (3) 


We apply Eq.(3) to compute the largest Lyapunov ex¬ 
ponent and FLI. We propose a new indicator called Ul¬ 
tra Fast Lyapunov Indicator (UFLI) with second order 
derivatives in order to detect chaoticity more rapidly and 
clearly as follows. The definition of UFLI is 
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where w n , Df{x n ), D 2 f(x n ) ( w n ) 2 ,w^ are a variational 
vetor, a Jacobian of f(x n ), a vector whose ith compo¬ 
nent consists of a product between *w n , Hessian ma¬ 
trix TL[fi(x n ))\ and w n where fi(x n ) is an ith compo¬ 
nent of f{x n ) and a orthogonal component of w n re¬ 
spectively. This proposal is different from the work by 
Dressier, Farmer [12] and Taylor [13] who introduce gen¬ 
eralized Lyapunov exponents using higher derivatives and 
the work by Barrio [6]. 


The formula (5) shows a variational equation consid¬ 
ering a second order derivative. We explain the ability of 
UFLI. Let us define u and inner product B™ by 

u n = D/„_iD/ n _ 2 • • • Df 0 w 0 , (7) 

Bi = Wft[/i(*n))]«n) (8) 
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Then, v n +i is expressed by 
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For example, at l = n, (n''_i DfjJ B®ei is considered 
as a variational vector evolved by Eq. (3) whose ini¬ 
tial vector is B®ei. Then, ^TIj=i -Bfe* expands 
in the direction of the Largest Lyapunov exponent. In 
the same way, higher order terms are considered as vari¬ 
ational vectors evolved by Eq. (3) if n is enough large. 
Therefore, ||«; n || grows drastically. The time evolution of 
UFLI changes clearly if an initial point belongs to chaos 
domain and grows slowly if the initial point belongs to 
a domain of KAM- or Resonant torus. Here, We ap¬ 
ply UFLI to the Froeschle map which is known to show 
Arnold diffusion and Chirikov diffusion [14-16], where 
the map is defined by 
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where I \, I 2 are action variables and Q \, 62 are action- 
angle variables correspondng to action variables respec¬ 
tively. Figures. 1, 2 and 3 show the time evolutions of 
UFLI, OFLI 2 , OFLI and FLI with the initial points A, B 
and C respectively. The floatl28 precision is used to cal¬ 
culate them. Three initial points A= (I \, 6 \, I 2 , 02 ) = 
(2.04,0,2.1,0), B= (1.8,0,1.2,0), C= (1.67,0,0.91,0) 
correspond to the chaotic domain, the KAM torus do¬ 
main and the resonant torus domain respectively in the 
Froeschle map with e = 0.6 [15]. We set the initial vari¬ 


ational vector as below. 


wi(0) 

= 0.001, 

w 2 (0) 

= 0.001, 

w 3 (0) 

- vs ~ 1 x 0.001. 

wj 4 (0) 

= 0.001, 

IM0)|| 

- x 0 001 


~ 0.0017. 


(13) 


According to Figs. 1, 2 and 3, our proposed UFLI per¬ 
forms much better compared to OFLI 2 , OFLI and FLI. 
Figures 4 and 5 show diagrams of UFLI and OFLI 2 for 
Froeschle map with e = 0.6 at n = 200 whose initial 
condition is 9 1 = #2 = 0. UFLI can show Arnold web, 
a structure consists of resonant lines more clearly than 
OFLI 2 . According to Ref. [16], this map behaves dif¬ 
ferently as a magnitude of e. It is known in Ref. [16] 
that , Arnold diffusion occurs when e < 0.9 and Chirikov 
diffusion occurs when e > 1.3. Here, we apply UFLI 
to detect a change between these diffusion regime. We 
compare variations of UFLI and OFLI 2 v.s. e. One thou¬ 
sand initial points are chosen near (Ji,^) = (7r/2,7r/2). 
Figure 6 shows ensemble average of UFLI(50) v.s. the 
parameter e change and Fig. 7 shows the counterpart of 
OFLI 2 (50), OFLI(50) and FLI(50) instead of UFLI(50). 
Figure 6 shows that UFLI loses smoothness in e > 0.9 
and distinguishes a transition between the two regimes 
(Arnold diffusion and Chirikov diffusion) of Froeschle 
map although OFLI 2 and the other existing detectors 
such as FLI cannot detect any transition in Fig. 7. We 
define Lyapunov time Tl and the macroscopic instabil¬ 
ity time T/. In Arnold diffusion, Tj ~ exp (Tl) while in 
Chirikov diffusion, Tj is order of polynomial for Tl [17]. 
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FIG. 1. Point A. A time evolution of UFLI, OFLI 2 , OFLI 
and FLI in a chaos domain. The common logarithm is used 
to calculate UFLI. OFLI 2 , OFLI and FLI. 
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FIG. 2. Point B. A time evolution of UFLI, OFLI 2 , OFLI 
and FLI in a KAM torus domain. The common logarithm is 
used to calculate UFLI, OFLI 2 , OFLI and FLI. 
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FIG. 3. Point C. A time evolution of UFLI, OFLI 2 , OFLI 
and FLI in a resonant torus domain. The common logarithm 
is used to calculate UFLI, OFLI 2 , OFLI and FLI. 
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FIG. 4. A digaram of UFLI for Froeschle map with e = 0.6 at 
n = 200. The common logarithm is used to calculate UFLI. 



FIG. 5. A diagram of OFLI 2 for Froeschle map with e = 0.6 at 
n = 200. The common logarithm is used to calculate OFLI 2 . 


Then, it is important to detect the point of e around 
which diffusion type changes. 


According to the result above, our proposed UFLI 
chaos detector is very powerful to detect chaoticity of 
systems with relatively small Lyapunov exponents more 
rapidly and clearly than FLI, OFLI and OFLI 2 . In addi¬ 
tion to this, UFLI can also detect a sharp change of the 
diffusion regime between Arnold diffusion and Chirikov 
diffusion although OFLI 2 and the other existing indica¬ 
tors cannot detect any transition. 

Log Fast Lyapunov Indicator Here, we investigate 
further to chaotic systems with zero Lyapunov exponent. 
In generally, a positive Lyapunov exponent shows a ex- 
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FIG. 6. The ensemble average of UFLI(50) v.s. e. The com¬ 
mon logarithm is used to calculate UFLI. 



FIG. 7. The ensemble average of OFLI 2 (50), OFLI(50) and 
FLI(50) v.s. e. The common logarithm is used to calculate 
OFLI 2 . 


istence of exponential growth of a variation between two 
close orbits. A positive value of Lyapunov exponent is 
used as an indicator of chaoticity. Even though the value 
of Lyapunov exponent is zero, behavior on torus and sub¬ 
exponential behavior with zero Lyapunov exponent are 
different. Thus, we propose another new indicator to 
distinguish them. In this section, Log Fast Lyapunov In¬ 
dicator (LFLI) 


LFLI(n) = log 


sup log 

0<i<n 


f \\Df(xi)wi\\ Y 

\ IKII ). 


(14) 


is proposed to characterize sub-exponential behaviors. 
Here, Wo,Wi,Df(xi) are an initial variational vetor, a 
variational vector at n = i and a Jacobian of f(xi) 
respectively. If infinite ergodic systems behave sub- 
exponentially [7], the time evolution of LFLI grows lin¬ 
early with slope smaller than unity. If systems have a 
positive Lyapunov exponent, the slope is unity. We apply 


LFLI to the Boole transformation and the S'-unimodal 
function in the following section. 

Boole transformation Here, the Boole transformation 
T : R —> R is defined by 


*£ 71+1 — T(x n ) — x n . (15) 

X n 

It is known that the Boole transformation is ergodic and 
preserves the Lebesgue measure [10]. The Boole trans¬ 
formation is an infinite ergodic system and the following 
equation is known to hold 

1 ” 

lim — Y f(T k x) = 0, (16) 

n—> oo 77 , * ^ 
k= 1 

a.e. x € R, v / € L 1 (/r), 

where p is the invariant measure for the probability 
preserving Boole transformation T [9]. By substitut¬ 
ing f{T k x) = log |T'(xfc)|, we know that the a value 
of Lyapunov exponent of Boole transformation is zero. 
However, it is known that the dynamical system behaves 
sub-exponentially [7]. Namely, its orbital expansion rate 
A grows A ~ exp(t 2 ). By using the LFLI, we can find a 
power index. To compare with the Boole transformation, 
we consider the following generalized Boole transforma¬ 
tions 


3*77+1 — Lfy , ta{Xn) — OiX n , (17) 

X n 

0 < a < 1,0 < P, 

which are known to have non-negative Lyapunov expo¬ 
nent [ 8 ] 


A = log ^1 + 2y/a(l - a) j . (18) 

We put /3 = a here for simplicity, because (3 doesn’t 
affect on the Lyapunov exponent A. Figure 8 shows en¬ 
semble averages of the time evolution of LFLI for the 
Boole transformation and the generalized Boole trans¬ 
formations with a = 0.99995 whose three hundred initial 
points are chosen near a point x = 11.7. Here, the ini¬ 
tial condition is ui( 0 ) = 0.00000000000000000001 and the 
float 128 precision is used to calculate LFLI. 


In Fig. 8 , the /(ln(n)) and g(ln(n)) are linear approx¬ 
imations of the ensemble averages of the Boole trans¬ 
formation and the generalized Boole transformations re¬ 
spectively. The slopes of /(ln(n)) and <?(ln(n)) are about 
0.436 and 0.957 respectively. These results indicate that 
our proposed LFLI is very powerful to find a power index 
for sub-exponential behavior. 
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FIG. 8. The ensemble average of the time evolution of LFLI of 
Boole transformation and generalized Boole transformation. 
The common logarithm is used to calculate LFLI. 


chosen near a point X = e _1 . Here, the initial condition 
is ut(0) = exp(—20) and the float 128 precision is used to 
calculate LFLI. 

Conclusion We propose two chaos indicators Ultra 
Fast Lyapunov Indicator (UFLI) and Log Fast Lya¬ 
punov Indicator (LFLI). It is found that UFLI can detect 
chaoticity more rapidly than OFLI 2 , OFLI and FLI and 
the only UFLI can detect a sharp change between Arnold 
diffusion and Chirikov diffusion regimes, that has not 
been detected by the existing methods such as OFLI 2 . 
LFLI can measure a power index of a sub-exponential 
system. In particular, LFLI firstly characterizes chaotic¬ 
ity of systems which have zero Lyapunov exponent which 
has been regarded as non-chaotic systems. Such detec¬ 
tors UFLI and LFLI proposed here are very promising 
to detect chaoticity of experimental data of intrinsically 
weakly chaotic systems. 
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FIG. 9. The time evolution of ensemble average of LFLI of 
the S-unmodal function. A slape of linear approximation for 
LFLI is about 0.486. The common logarithm is used to cal¬ 
culate LFLI. 


S-unimodal function S'-unimodal function [11] is de¬ 
fined by 

1 c V 2 

X n +1 = Ts(X n ) = 1 +3X 2 , Xn € [-1,1]. (19) 


S'-unimodal function has an infinite measure whose den¬ 
sity function is 


P(X) 


1 

(l + x)VT=^*' 


( 20 ) 


According to Ref. [11], orbit expansion rate of S- 
unimodal function is A ~ exp(n5). Figure 9 shows the 
ensemble average of the time evolution of LFLI for the S- 
unmodal function whose one thousand initial points are 
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